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ABSTRACT 


The atmospheric drag on a satellite is conveniently determined from the 
rate of change of the orbital period dP/dt. This procedure can be extended 
to the near-free-molecule-flow regime, in which the drag coefficient of a 
sphere is expressed by C D = C Dfm - F(S b ,S 00 )/K, where C Dfm is the free- 
molecule-flow coefficient, and F(S b , S^) is a function of the molecular speed 
ratios for gas diffusely emitted from the body surface and for the free stream. 
The free-stream Knudsen number K can be expressed as 1 /K = BDp, where 
D is the sphere diameter, p is the atmospheric density, and B is a constant. 
In terms of the usual orbital elements and the density at perigee pp, 

dP/dt= - 3a (s)( £ l C DfmPp- f 2 BDF 4) • 


where 


TT 

f = f /_£_{ (1 + e cos E) 3/2 
1 J \ p P/ (l-ecosE) 1/2 


dE 


From the measured dP/dt for two or more satellites with nearly identical 
A/m and orbits, it would be possible to find the atmospheric density at low 
altitude and a value of F for comparison with theory. 

For a practical experiment, particular orbit characteristics and satellite 
parameters must be selected to match available launch-vehicle performance. 
Two concrete examples are discussed. The first assumes the performance 
of the Scout vehicle. This example shows that a reasonable experiment 
seems possible even for this modest vehicle. The second example adopts 
parameters that may be typical of a secondary mission on a Saturn class 
vehicle. 
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DESIGN OF A SATELLITE EXPERIMENT 
FOR ATMOSPHERIC DENSITY AND NEAR-FREE-MOLECULE-FLOW 

AERODYNAMICS 


Louisa S. Lam, Geraldine M. Mendes, and Charles A. Lundquist 

1. INTRODUCTION 

Accurate tracking of artificial satellites has provided a powerful method 
for studying the several phenomena that result in forces acting on an orbiting 
body. Meticulous analyses of satellite positions have produced a detailed 
representation of the earth's gravitational potential (Gaposchkin, 1966a, 1967; 
Kozai, 1964, 1966) and models for the atmospheric density and its variation 
(Jacchia, 1965). Most of the atmospheric information obtained this way per- 
tains to altitudes above 200 km. It is natural to inquire how established pro- 
cedures can be extended to study the atmosphere and aerodynamic phenomena 
at lower altitudes for which departures from free-molecule-flow aerodynamics 
are significant for bodies of practical dimension. These altitudes are below 
those where satellites usually persist for appreciable lifetimes. 

Within the possibilities of present technology, reasonable experiments 
can be designed that will yield both a measurement of atmospheric density 
and experimental data on the drag coefficient of a sphere in near-free- 
molecule flow. This dual objective is natural, even inevitable, since the 
aerodynamic forces determined from resulting satellite accelerations, of 
course, depend upon both the atmospheric density and the details of the inter- 
action of the satellite with the atmosphere. 


This work was supported in part by NsG 87-60 from the National Aeronautics 
and Space Administration. 
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The suggested experiments require that a single rocket vehicle launch 
an ensemble of two or more spherical satellites into essentially identical 
orbits. These bodies would be similar in most of their characteristics, but 
would differ in those that allow a separation and an identification of the factors 
that enter into the aerodynamic drag. 

In arriving at a practical experiment, we must recognize the satellite - 
tracking accuracy currently attainable. The accuracy obtained with the exist- 
ing network of Baker-Nunn cameras has been assumed. 
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2. AERODYNAMIC-DRAG REPRESENTATION 


To design the experiment, we must adopt a representation of the drag 
coefficient for near-free-molecule flow. The choice is somewhat subjective, 
since several formulas have been suggested by different authors. Fortunate- 
ly, some of the earlier differences between theory and laboratory results 
have been resolved in later papers (Maslach, Willis, Tang, and Ko, 1964). 
Nevertheless, further data of the sort obtainable from a satellite experiment 
would be particularly valuable to compare with theory for parameter values, 
such as large Knudsen numbers, difficult to attain in the laboratory. 

Most of the theoretical expressions for the near-fr ee-molecule-flow 
drag coefficient of a sphere have the form (Maslach et al. , 1964) 


C D= C D£m-I r ' S b' S J • (1) 

where is the value of the drag coefficient for fr ee-molecule flow, 

K is the free-stream Knudsen number (free-stream mean free path/sphere 

diameter), is the molecular speed ratio of the gas diffusely emitted from 

the surface, and S is the molecular speed ratio for the free stream. For 

oo 

this discussion, we can represent the last three quantities by the following 
equations : 


K = 


w 


n/T IT N n (T p D 


( 2 ) 


s 


b 


u 

V 2RT b /w 


( 3 ) 


3 


and 


S 

00 


U 

x /2 RT 7w 
\ 00 


(4) 


where D is the sphere diameter, p is the atmospheric density, w is the 

mean molecular weight of the air molecules, N„ is Avogadro's number 
23 u 

(Nq = 6. 02 X 10 atoms / g mole), <r is the effective collision diameter of a 

molecule of air, is the surface temperature of the orbiting body, 

is the temperature of the atmosphere, U is the satellite velocity 
relative to the atmosphere, and R is the gas constant per mole 

7 

(R = 8.3170 X 10 erg/mole °K). In this same notation, for diffuse reflec- 
tion (Schaaf and Chambre, 1958), 


C Dfm 


2 + 


1.18 


(5) 


Different forms of the function F(S, , S ) have been derived by various 

b oo 

authors using different procedures (Maslach et al . , 1964): 


(Baker and Charwat) F = 0. 24 + 1 . 06 , (6a) 

0. 165 S + 1.44 - 1. 13/S 

(Willis) F = , (6b) 

00 

0. 33 S - 0. 12 

(Rose) F = g . (6c) 

00 


In their review of this topic, Maslach et al . conclude that the last two expres- 
sions give the best agreement with laboratory data. 
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Using equation (Z) to express 1/K in terms of p, w, and D, we get a 
representation of the aerodynamic drag per unit mass of the satellite: 


drag J_ c A_ 2 

unit mass Z Dm 


p = 2 


C Dfm - BF(S b’ S J D P 


A. tt 2 
— p U 
m r 


( 7 ) 


where A is the cross-section area of the sphere, m is the mass of the sphere, 
and 

\[2 IT 0- 2 N n 

B = . 

w 

For the purpose of the experiment under discussion, we can consider 

several of the quantities in this representation to be known. The velocity of 

the satellite can be found from the orbital elements obtained by tracking the 

satellite. The surface temperature of the sphere can be calculated using 

proved techniques (e. g., Heller, 1961). The mean molecular weight and the 

effective collision diameter of the air molecules vary quite slowly with 

altitude in the interval of interest, and appropriate values from the literature 

- 1 0 

can be used ( cr - 3. 65 X 10 m). The diameter and mass of the sphere 

are parameters to be adjusted in the design of the experiment, but will be 
known quantities for an actual body. 

If equation (5) is accepted as correct, only the atmospheric density p 
and the values of F(S^, ) remain to be determined from appropriate analy- 

sis of tracking data. The value of p is of obvious geophysical interest. The 
experimental values of F(S^, S ) are valuable for comparison with theoretical 
values derived by various procedures (equations (6)). 

To estimate the magnitude of the second term of equation (7) relative to 
the first term, we can use the theoretical expressions for F with values of 
p from a standard atmosphere (Jacchia, 1965). Values for equation (6b) 
are shown in Table 1. 
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Table 1. Drag coefficient for 1-m-diameter sphere 


Altitude (km) 

K 

BF < s b- s J D 0 

C Dfm 

C D 

120 

3. 1 

0.084 

2. 062 

1.978 

130 

9.4 

0. 035 

2. 063 

2. 028 

140 

19. 2 

0. 019 

2. 063 

2. 044 

150 

32. 6 

0. 012 

2. 064 

2. 051 

160 

50. 2 

0. 009 

2. 064 

2. 056 

170 

72. 2 

0. 006 

2. 065 

2. 059 

The following values are used in the table: 


atmospheric model from Jacchia (1965) with exos 
temperature = 1300°K, 

F is given by equation (6b), 

T b = 300° K, 

U =8.2 km/sec. 

pheric 
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3. ORBIT CHANGES RESULTING FROM DRAG 


The effect of atmospheric drag on the anomalistic period P of a satellite 
has been conveniently expressed by Sterne (1958). If the rotation of the 
atmosphere is neglected, the rate of change of the period is 


AP _ dP 
P " dt 


1 A. 

2 m 




_£_ il_L e cos E)^ ^ 
Pp (1 - e cos 


dE 


(8) 


where Pp is the density at perigee, a the semimajor axis, e the eccentricity, 
and E the eccentric anomaly. As explained by Jacchia (1963), this equation 
is the usual starting place for determinations of atmospheric density because 
the period and its changes are rather directly obtained from tracking data. 


To simplify notation. 


it is convenient to define quantities f^ and f^ by 


and 



IT 



(1 + e cos E) 31 * dE 
P P (1 - e cos E) ' C 



3/2 

( 1 + e cos E) 

(1 - e cos E) 1 /2 


dE 


(9) 


( 10 ) 
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Tables of are given by Jacchia and Slowey (1963a) and similar tables of f? 
are easily constructed. Using this notation and the adopted expression for 
Cp in equation (8) gives (Lundquist, 1967) 

¥=- 3 af.A. r p 0 + 3 af, B — D F p* . (11) 

dt 1 m Dfm r P 2 m 

In this equation, B has been treated as a constant because it varies much 
more slowly with altitude than does p. 
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4. CONCEPT OF EXPERIMENT 


To investigate near -fr ee-molecule-flow aerodynamics and atmospheric 
density at low altitudes, we suggest two spherical satellites having identical 
A/m and having surface preparations that ensure equal surface temperatures. 
These spheres should be launched by the same vehicle into initial orbits as 
nearly identical as practical. The diameter of one sphere, D(l), should be as 
small as possible and still allow reliable photographic tracking with the 
Baker-Nunn cameras. The diameter of the second, D(2), should be as large 
as can be contained within the available launch vehicle. Using simple proce- 
dures discussed below, values for the rate of change of the periods P(l) and 
P(2) can be obtained from tracking data. For two satellites under these con- 
ditions, equation (11) gives two linear simultaneous equations (12) and (13) 

2 

for unknowns C-p£ m Pp an ^ BFpp . The coefficients multiplying these unknowns 
involve only quantities known a priori or determined from tracking data. 

The quantities f^ and f^ are functions of the atmospheric scale height near 
perigee, as well as of the orbital elements, but the value of the scale height 
can be estimated with sufficient accuracy for computing f^ and f^r 


- 3 


, A 
af 7 — 
1 m 


C Dfm P P + 


3 af D( 1 ) — 
2 m 


BF Pl 


= dP(l) 
dt 


( 12 ) 


and 


- 3 


, A 
af, — 
1 m 


C Dfm P P + 


3 af D(2) — 
2 m 


BFp?, = 

p P dt 


(13) 
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Under the assumption that the orbits are identical, these have the simple 
solutions 


C Dfm P P 


1 

( 3 af i [D(2)-D(l)] 


D(J) 


dP(2) dP(l) 

~d^" ~dt~ 


( 14 ) 


and 


BFPp 


1 

dP(2) 

dP(i) 

(3 af 2 [ D (2) - D(l)] 

dt 

dt 


(15) 


If the perigee altitude were above the near-free-molecule-flow region 
for the adopted satellite diameters, the periods would have the same rate of 
change and equation (14) would reduce to the usual expression applied to 
determine Cp^^Pp* Thus equation (14) is only the familiar expression with 
a simple correction for the case of near - fr ee-molecule flow. If a value of 
^Dfm * S acce P te< ^ f rom equation (5) or from some other source (Cook, 1965), 
equation (14) yields a value for p . Using this value for p with equation 
(15) gives a value for BF, which is a measure of the departure from free- 
molecule-flow aerodynamics. If B is considered known, a value of F follows 
immediately. 

The accuracy with which F is obtained depends critically upon the selection 
of satellite and orbit parameters such that the rates of change of the periods 
differ significantly. Variation of pp during the interval of observation can 
also limit the accuracy obtained. 

Equations (14) and (15) are intended only to illustrate the principal fea- 
tures of the analysis. For example, as the satellites separate, the assump- 
tion of identical orbits begins to break down; however, first-order correc- 
tions for this can be introduced. The effects of an initial separation 
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velocity between the spheres and the rotation of the atmosphere can likewise 
be incorporated. 

Introducing into the experiment a third sphere, again with the same 
A/m but with a still different diameter, D(3), would provide a useful check 
on the results. A third equation of the form (12) is not independent, but it 
must be consistent with (12) and (13). A demonstration of the required consis- 
tency would lend substantial credence to the solution. This is particularly 
true because equations (12) and (13) will always have a solution, but this solu- 
tion is by no means guaranteed to be meaningful, as it assumes equation (1) 
is correct. Given three or more equations of the form of (12) with reasonable 
consistency, a least- squares solution can be obtained replacing equations 
(14) and (15). 

On the other hand, there is some reason to believe that a term proportional 
2 

to (D pp) should be added to equation (1) (Sherman, Willis, and Maslach, 1964). 
If this should be necessary, equation (11) would have a corresponding ad- 

3 

ditional term involving pp, f^, and a coefficient characterizing the importance 
of this added term. Data on three spheres of different size then would be 
required to investigate the departure from f ree-molecule-flow aerodynamics, 
but otherwise the same procedures could be followed. 
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5. EXTENDED EXPERIMENT 


The simple experiment sketched in the previous section can be extended 
to alleviate the weaknesses that result because the orbits become different as 
the spheres separate. This is accomplished by adjusting slightly the masses 
and hence the A/m ratios to just balance the expected drag differences caused 
initially by departure from fr ee-molecule-flow aerodynamics. For the two- 
sphere case, this requires that (A/m)(l) and (A/m)(2) be chosen so that 


dP(l ) _ dP(2) 
dt dt 


( 16 ) 


This implies that 

> f , C Dfm - 

—(2) 'hDfm-VX^P ' ‘ ’ 

nr 7 

Because BF and Pp are imperfectly known, equation (17) cannot be 
exactly satisfied; however, the best estimates of these can be used to achieve 
a first approximation. Also, as the eccentricities of the orbits decrease, the 
values of f ^ and f^ change, which may generate a departure from equation 
(17). 


Corresponding to equations (14) and (15) for this generalized case are the 
equations 

c = 1 

Dfm pp [3 afj ~(1 ) “(2)J [D(2) - D(1 )] 

( 18 ) 




D(2) —(2) d ^ 2) 
m dt 


12 



and 


BF Pp = 


1 


L 3 af 2 ^(l)-(2)][B(2) -D(l)] 




- —( 2 ) 

m 


dP(l) 

dt 


(19) 


Perhaps it is worth noting that these equations are meaningful whether or 
not equation (17) is satisfied. 


In the following discussions and examples, the point of view of the pre- 
vious section will be retained because the spatial separation of the spheres 
is an easily pictured measure of the observability of the desired quantities. 
However, an actual execution of the experiment would probably employ the 
more general approach of this section. 
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6 . TEMPORAL VARIATION OF THE ATMOSPHERE 


Changes of atmospheric density as a function of time, are well known 
above 200 km, but less is known of their magnitude at lower altitudes 
(Jacobs, 1967). Measurements of changes that may occur are as interesting 
as determinations of the density. An investigation of temporal phenomena 
would be facilitated by the inclusion of a very heavy sphere in the ensemble. 
This sphere could have the same size and surface preparation as one of 
the other spheres, but it should have as great a mass as practical. 

The high-density body will remain in orbit longer than the lighter bodies 
in the aerodynamics experiment. The knowledge of the drag coefficient 
generated in the initial phases of the experiment is directly applicable to the 
later history of the heavy body. An analysis of the orbit of the long-lived 
satellite in the usual way is useful in the measurement of the variation of 
atmospheric density as a function of time and of solar and geophysical 
events . 


14 



7. SURFACE PREPARATION 


In the adopted expressions for the drag coefficient, the properties of the 
surface enter through the ratio of the orbital speed to the speed with which 
gas molecules are diffusely emitted from the surface. The latter speed de- 
pends upon phenomena at the surface of the satellite. The expression for 
given in equation (3) is an approximation for more involved theories that 
could be invoked. For example, the relation between the surface temperature 
and the velocity of emitted molecules is often characterized by an accom- 
modation coefficient (Wachman, 1962; Springer and Tsai, 1964). 

For the initial experiment considered here, the most important require- 
ments are that the spheres be carefully prepared with identical surface 
characteristics with respect to velocity of emitted molecules, and that the 
surface temperatures be equal. It is hoped that these conditions are not 
contradictory. These requirements have been assumed in the description 
of the experiment. Alternatively, the theory of the experiment might be 
slightly modified to include a correction for unavoidable differences in 
surface properties. 

The usefulness of the experiment can be enhanced if the surface prepara- 
tion is one that has been studied by laboratories investigating the relevant 
surface phenomena. Information from such studies will be valuable in the 
selection of an appropriate surface material and preparation. 

If vehicle performance would allow additional satellites, an experiment 
with spheres differing only in their surface characteristics could be 
envisioned, and the data obtained could be compared with predictions from 
theory. There is little reason, however, to perform such an experiment 
at very low-perigee altitudes, unless the objective were a study of the depend- 
ence of F on S^. The latter would presuppose an adequate evaluation of 
S b and the dependence of C^ m upon it. 
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8. REPRESENTATION OF THE GEOPOTENTIAL 


The determination of accurate orbits for the satellites in the suggested 
experiment depends upon an adequate representation of the geopotential. 
Fortunately, the philosophy of the experiment minimizes the demands on the 
geopotential representation. Because the perigee of the suggested satellites 
will be quite low, and if the eccentricity is not too great, the effect of high- 
order harmonics in the gravitational field may be greater than that experienced 
by satellites in more distant orbits. Thus the adopted standard representation 
of the geopotential may not account for all detectable gravitational perturbations 
on an individual satellite. This may introduce an increased "noise level" in 
the orbit-determination process. 

If this should become troublesome, the same dense sphere used for 
temporal atmospheric variation also may relieve this problem. If the 
additional sphere is denoted by 4, then from equation (8), 


dP(4) = m(l) dP( 1 ) 
dt m(4) dt 


( 20 ) 


If, further, m(4) is significantly larger than m(l), perhaps by a factor of 5 or 
10, then equation (20) provides accurate values of dP(4)/dt from the larger, 
more easily measured values of dP(l)/dt. 

A similar situation prevails for the drag-induced changes in the other 
orbital elements of sphere 4. These can be expressed with relatively greater 
accuracy from the drag measured by sphere 1. But the drag on sphere 4 is 
already relatively small because A/m is small. Hence, the uncertainty 
associated with drag on sphere 4 is something like m(l)/m(4) less than 
that on sphere 1 . 
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Drag corrections derived from sphere 1 can be applied to the orbit 
of sphere 4, and this orbit examined for the adequacy of the adopted geo- 
potential. If the residuals between observation and calculation approach 
the accuracy of the observations, then the adequacy of the geopotential is 
strongly verified. In this case, the analysis of atmospheric effects on the 
light satellites can proceed with assurance that the ’’noise 11 induced by uncer- 
tainty in the geopotential is small. Asking that the residuals approach the 
observational accuracy is probably too stringent. Presumably, the geopoten- 
tial is sufficient if the residuals of the corrected orbit for sphere 4 are small 
compared to the uncertainty in the determination of the separations induced 
by differences in for the high-drag spheres. 

If the residuals for the corrected orbit of sphere 4 are greater than 
expected, the orbit can perhaps provide information for an improved deter- 
mination of the geopotential. In an imprecise sense, the corrected orbit is 
a low-perigee, drag-free orbit. 
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9. OBSERVATION AND DATA PROCESSING 


The Baker-Nunn cameras and associated timing systems produce 
fundamental data in the form of station-to-satellite directions accurate to 
2 to 3 arcsec and associated observation times accurate to 1 msec. The 
station positions are known globally to an accuracy of 10 to 20 m (Kohnlein, 
1966, 1967). The observations and station positions are used in a differential 
orbit improvement (DOI) program to obtain an orbit for the satellite observed 
(Veis and Moore, I960; Gaposchkin, 1966b). This program uses a standard 
set of station coordinates and coefficients for the spherical harmonic repre- 
sentation of the earth's gravitational field (those currently used are from the 
Smithsonian Institution 1966 Standard Earth (Eundquist and Veis, 1966)). 

It also takes into account the lunar perturbations and various minor correc- 
tions. For a satellite on which drag is insignificant, the mean residual 
between observation and calculation is typically around 4 arcsec, which 
approaches the accuracy of the fundamental observations. 


When used on a satellite with appreciable drag, the DOI program gives 
a series representation for the orbital elements, including mean anomaly, 
in which the most significant remaining effect is that due to atmospheric 
drag. For the mean anomaly M this series is usually expressed in the form 

M = M q + Mjt + M 2 t 2 + . . . , (21) 

where are constants determined by the DOI program, and t is the time 
interval from some initial time. This is the form in which SAO publishes 
satellite orbital elements. In principle, dP/dt due to drag may be obtained 
directly from 
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In practice, a more refined but conceptually equivalent procedure is actually 
used (Jacchia and Slowey, 1963b). The values obtained are available for use 
with equations (14) and (15), or with (18) and (19). 
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10. EXPERIMENT REQUIREMENTS 


For a practical experiment, particular orbit characteristics and satel- 
lite parameters must be selected to match available launch-vehicle per- 
formance. Two concrete examples will be discussed. The first assumes 
the performance of the Scout vehicle. This example shows that a reason- 
able experiment seems possible even for this modest vehicle. The second 
example adopts parameters that may be typical of a secondary mission on a 
Saturn class vehicle. 

In the Scout example, only two spheres are considered, recognizing the 
vehicle limitations. In the second example, three spheres are suggested for 
the aerodynamic experiment and a fourth sphere is included to monitor time 
variations of the atmosphere and to verify the geopotential representation. 
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11. EXAMPLE 1, SCOUT LAUNCH VEHICLE 


For this example the performance of the Scout vehicle is adopted (National 
Aeronautics and Space Administration, 1966). Within necessary constraints, 
the objective is to define an optimum experiment. The first constraint, im- 
posed by the philosophy of the experiment, requires that the two spheres have 
equal area-to-mass ratios, so that 


D(1) Z = D(2) 2 

m(l) m(2) 


(23) 


For an orbit yet to be specified, the sum of the two masses must correspond 
to m, the Scout payload capability less an allowance for support structure. 


m(l) + m(2) = m 


(24) 


The diameter D(2) of the larger sphere will be fixed at the maximum 
diameter that can be contained within the Scout payload area. This is about 
76 cm (National Aeronautics and Space Administration, 1966): 

D(2) = 76 cm . (25) 


To ensure that the small sphere can be tracked with the Baker-Nunn 
camera, a minimum diameter must be imposed. This minimum depends 
upon the maximum distance to the satellite. This condition can be expressed 
as a function of the apogee distance r^ and the radius of the earth r : 

D(l) = 10" 7 (r A - r E ) . (26) 

A linear relationship between diameter and range follows from the usual 
expressions for illumination from a satellite as a function of diameter and 

_7 

radius (Veis, 1965; Zirker, Whipple, and Davis, 1956). The factor 10 is 

a threshold based on experience with the Baker-Nunn camera. A more con- 

-7 

servative value would be 1. 4 X 10 
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Combining (23) through (26) gives m(l), m(2), D(l), and D(2) as functions 
of the parameters related to the vehicle performance, i.e., m, Tp, and r^. For 
an eastward launch from Wallops Island, Scout performance data are avail- 
able as a family of curves of apogee versus payload for a sequence of perigee 
values. Hence, for various values of Tp and r^, a value of m is obtained, and 
consequently values of m(l), m(2), D(l), and D(2) can be calculated. 

A set of these can be used with a theoretical value of F and a nominal 
value of pp in equation (15) to predict a value for [dP(2)/dt] - [dP(l)/dt] . 
Curves obtained this way are shown in Figure 1. The Jacchia (19&5) model 
atmosphere with an exospheric temperature equal to 1300°K has been used for 
these and subsequent examples. 

The difference in the rate of change of the periods is a measure of the 
detectability of the departure from free-molecule-flow aerodynamics. A 
more easily interpreted measure is the separation of the two satellites 
after some time interval. This can be estimated starting with equations (12) 
and (13) and using the theoretical values of the required parameters. 

Assuming for this estimate that there is negligible initial separation velocity, 
the expected separation in mean anomaly can be estimated by the equation 

AM= =2l FB i p P ^[D(2)- D ( 1} ]t 2 , (27) 

where AM is in revolutions, and P is the mean period of the two satellites 
for the interval t since launch. 

If the atmospheric density at a particular altitude is an objective of the 

experiment, this altitude determines Tp. In this case, the corresponding 

curve in Figure 1 shows the range of obtainable r^ and [dP(2)/dt] - [dP(l)/dt], 

still assuming Scout performance. Equation (27) can be used to estimate the 

separation of the satellites as a function of time. In Figure 2, the factor 
2 

multiplying t on the right side of equation (27) is plotted for various values 
of r p and r^. The initial value of P has been used for P. Recognizing this con- 
servative approximation, the separation at time t after launch can be found. 
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A rough lower limit on the useful separation is set by the tracking 
accuracy; the separation should be many times greater than 2 arcsec as 
seen from the earth's surface. A hazy upper limit corresponds to a separa- 
tion so great that the satellites are no longer in nearly identical orbits. The 
film dimension of the Baker-Nunn camera (30° topocentric) provides an alter- 
native limit that might be adopted if both satellites were to appear on the same 
film; this is hardly necessary, however. 

Another sort of limit is imposed by the time interval available for 
observation. In Figure 1, satellite orbital lifetimes have been noted at 
intervals along the curves. The lifetimes have been estimated using the 
procedure described by Jacchia and Slowey (1963a). Probably, observations 
should be made over at least a month for a statistically sound orbit deter- 
mination. For practical reasons, it may be reasonable to limit the duration 
of the experiment to a few months. Surely a duration greater than a year is 
dubious. 

A specific example will illustrate the open possibilities. Suppose that 
the atmospheric density at 150 km is an objective and that 3 months is 
adopted as a reasonable duration for the experiment. From Figure 1, an 
apogee altitude of about 2400 km is indicated for these conditions. Also, 
Figure 1 gives [dP(2)/dt] - [dP(l)/dt] = 2.3 x 10 for the initial relative rate 
of change of the period. From Figure 2, the spheres would separate about 
0. 1 rad during the first month. The separation would be some 0.4 rad 
after 2 months and would be something like 1.0 rad before the satellite 
demise after 3 months. The satellite diameters, masses, and orbital param- 
eters for this case are given in Table 2. 

This example illustrates an important characteristic of the orbits satis- 
fying the imposed conditions — namely, they are initially quite eccentric. 

This is fortunate for it facilitates reliable tracking. 
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Table 2. 


Orbital parameters at apogee 2400 km 


rp - rg = 1 50 km 
r^ - Tp = 2400 km 
a = 7653 km 
e = 0. 147 
m = 1 19. 3 kg 
D (1 ) = 24 cm 
D (2) = 76 cm 
m( 1 ) = 10.8 kg 
m(2) = 108. 5 kg 

2 -1 

A/m = 0.0418 cm g 
-3 

Satellite density (1) = 1.49 g cm 

-3 

Satellite density (2) = 0.47 g cm 
Lifetime = 90 days 


An alternative criterion for selecting a set of satellite and orbit param- 
eters uses equation (27) with the t replaced by the lifetime of the satellite. 
The resulting value AM is a figure of merit akin to the total separation of 
the satellites during their life. Using the same Scout performance data. 
Figure 3 shows this AM as a function of apogee for the perigee altitude of 
150 km. In some sense, the maximum of such a curve represents an opti- 
mum situation. The characteristics of this case are shown in Table 3. 
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Figure 3. Separation as a function of apogee, for example of Scout- 
launched satellites. 
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Table 3. 


Orbital parameters at apogee 4000 km 



a 

e 

m 


D(l) 

D(2) 

m(l) 

m(2) 

A/m 

Satellite density (1) 
Satellite density (2) 
Lifetime 


= 150 km 

= 4000 km 

= 8453 km 

= 0. 228 

= 99.5 kg 

= 40 cm 

= 76 cm 

= 21.9 kg 

= 77. 6 kg 

= 0. 0573 cm 2 g" 1 
_ 3 

= 0. 654 g cm 
= 0. 338 g cm" 3 
= 131.9 days 
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12. EXAMPLE 2 , SATURN SECONDARY EXPERIMENT 


For a secondary experiment on a Saturn or another large launch vehicle, 
a set of nominal parameters is given in Table 4. The smallest sphere diam- 
eter again corresponds to the lower limit for reliable tracking, the largest 
diameter is dictated by available space within the vehicle, and the third 
sphere has an intermediate diameter. 


Table 4. Nominal parameters for Saturn secondary experiment* 


Satellite 

Diameter (m) 

Mass (kg) 

_ 3 

Density (g cm ) 

A/m (cm^ g *) 

i 

0 . 2 

15 

3. 58 

0. 0209 

2 

0 . 6 

135 

1.19 

0 . 0209 

3 

1.0 

375 

0. 72 

0 . 0209 

4 

0.6 

612 

5.4 

0. 0046 

* 

perigee 

= 16 0 km. 




apogee 

= 2000 km. 




a = 7458 km. 




T, = 300° K. 
b 





Perigee altitude is selected as an interesting region of the atmosphere. 
The nominal apogee altitude then corresponds to about an orbital lifetime of 6 
months. The total mass of the ensemble corresponds to an estimate of vehicle 
payload that might be available for a secondary experiment (2500 pounds). 

In Table 4 , spheres 1, 2, and 3 have the same A/m. Sphere 4 has the 
same size as sphere 2, but a much smaller A/m. This last sphere has two 
purposes: detection of temporal changes of the atmosphere over an extended 
time interval, and examination of the adequacy of the geopotential represen- 
tation. 
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For the first three spheres, the relative rates of change in the period 
given by equation (12) are shown in Figure 4. Curves are given for a family 
of perigee and apogee values in the neighborhood of the nominal orbit in 
Table 4. 

For the cases in Figure 4, Figure 5 gives the corresponding coefficients 

2 

of t in equation (27). In Figures 6 and 7, the separation of the spheres as 
a function of time is illustrated for an apogee of 2000 km and various perigees. 
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(a) 


(b) 




k" r E (km) 


Figure 4. (a) - versus apogee for the example of Saturn- 

launched satellites. 

(b) 111!!- versus apogee for the example of Saturn- 

launched satellites. D(l) = 20 cm; D(2) = 60 cm; 

D(3) = 1 00 cm. 
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Figure 5. Curves for the example of Saturn- luanched satellites. 

versus apogee for: (a) sphere 2 from sphere 1, and (b) sphere 
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3 from sphere 1 . 
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Figure 6 
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Angular separation of sphere 2 from sphere 1 at apogee 
2000 km versus time, for Saturn-launched satellites. 
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F igure 7 
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Angular separation of sphere 3 from sphere 1 at apogee 
2000 km versus time, for Saturn-launched satellites. 




13. CONCLUSIONS 


In summary, an interesting class of experiments can be defined that 
achieve a dual objective. They measure the atmospheric density in the 
altitude range sometimes referred to as the "ignorosphere, " and they 
provide an experimental check on the theory of near-free-molecule-flow 
aerodynamics. 

The means for collecting the observational data are completely in 
existence in the SAO network of Baker-Nunn cameras. The data-analysis 
techniques are closely analogous to proved procedures used already to 
obtain much of the present knowledge of the atmosphere above 200 km. 

The launch-vehicle requirements are minimal for a satellite experiment. 
The Scout vehicle can accommodate a useful experiment. The complexity 
of the payload is also minimal, which makes the experiment attractive for 
a secondary payload on a launch vehicle. An ensemble of spherical satellites 
can be specified for more or less ambitious investigations, depending on the 
available vehicle capacity. 
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NOTICE 


This series of Special Reports was instituted under the supervision 
of Dr. F. L». Whipple, Director of the Astrophysical Observatory of the 
Smithsonian Institution, shortly after the launching of the first artificial 
earth satellite on October 4, 1957. Contributions come from the Staff 
of the Observatory. 

First is sued to ensure the immediate dissemination of data for satel- 
lite tracking, the reports have continued to provide a rapid distribution 
of catalogs of satellite observations, orbital information, and prelimi- 
nary results of data analyses prior to formal publication in the appro- 
priate journals. The Reports are also used extensively for the rapid 
publication of preliminary or special results in other fields of astro- 
physics. 

The Reports are regularly distributed to all institutions partici- 
pating in the U. S. space research program and to individual scientists 
who request them from the Publications Division, Distribution Section, 
Smithsonian Astrophysical Observatory, Cambridge, Massachusetts 
02138. 



